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Abstract 



O 

We have examined the dynamical behavior of the kink solutions of the one- 
J^ ■ dimensional sine-Gordon equation in the presence of a spatially periodic para- 

metric perturbation. Our study clarifies and extends the currently available 
knowledge on this and related nonlinear problems in four directions. First, 
/\ ' we present the results of a numerical simulation program which are not com- 

d ' patible with the existence of a radiative threshold, predicted by earlier calcu- 

lations. Second, we carry out a perturbative calculation which helps interpret 
those previous predictions, enabling us to understand in depth our numerical 
results. Third, we apply the collective coordinate formalism to this system 
and demonstrate numerically that it accurately reproduces the observed kink 
dynamics. Fourth, we report on a novel occurrence of length scale competi- 
tion in this system and show how it can be understood by means of linear 
stability analysis. Finally, we conclude by summarizing the general physical 
framework that arises from our study. 
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I. INTRODUCTION. 

Technological progress has made possible the fabrication of highly ordered materials and 
structures for a very large number of applications. In parallel to those advances, it has also 
been realized that the special properties required for many purposes necessitate inhomo- 
geneous systems. Here inhomogeneity may mean spatial modulations, quasiperiodicity, or 
disorder of several kinds. In addition, there are other situations in which inhomogeneity 
is undesirable but unavoidable. In either case, the study of disordered systems acquires 
fundamental importance. This is even more so when the physical system in which disorder 
or inhomogenity is to be studied is described by a nonlinear model. Whereas the role of 
disorder in linear problems is at least partially understood, much less is known about nonlin- 
ear disordered systems. In fact, even from a mathematical viewpoint, the understanding of 
these models, often related to stochastic partial differential equations (PDE), is very limited. 
Consequently, a great deal of research has been devoted to this topic [|^-|^ . 

A major part of the work done so far regarding nonlinear disordered systems has been 
concerned with some particular examples that are amenable to analytical treatment while 
capturing some essential physics. The sine-Gordon (sG) (actually, the whole family of non- 
linear Klein-Gordon equations, including, e.g., the 0^, double- and quadratic- sine-Gordon 
equations) and nonlinear Schrodinger (NLS) equations are often chosen as very suitable 
"canonical" examples. This is due to the fact that the basic mathematical structure under- 
lying them is well known and therefore provides a good starting point for theoretical work. 
This reason would not be sufficient if these models were not also related to a large number 
of phenomena that occur in quasi-one-dimensional physical systems, as is in fact the case. 
In the context of these two models, disorder is introduced through suitably chosen pertur- 
bation terms (see for an extensive list of physically relevant perturbations). This is the 
usual procedure by which inhomogeneity of any kind is studied: The equation describing the 
problem is established, the terms relevant to the considered physical situation are identified, 
and a perturbation to those terms is introduced, representing the desired kind of disorder. 
Our viewpoint in this work is more generic: Although the system we deal with is indeed 
related to a number of applications, our aim is that we will be able to gain insight into un- 
derlying mechanisms of the phenomenology of nonlinear disordered systems. Therefore, we 
introduce a simple periodic perturbation which will allow us to study very interesting and 
general phenomena, such as length scale competition, and will provide information relevant 
to the more complicated processes occurring in random media (the periodic potential can 
be interpreted as a "color" of a general noisy one). The knowledge obtained will also be 
useful to tackle other problems where detailed studies including analytical treatment are not 
possible. 

In this paper we study the behavior of one-dimensional (ID) sG kinks when perturbed 
parametrically by a spatially periodic potential. Initially, this was motivated by our related 
research (from the above general point of view) on the sG |^|5[ and NLS ||^J^ models. As a 
preliminary step to the investigation of sG breather dynamics Q on these kind of potentials, 
it is natural to first seek a good understanding of kink dynamics. Therefore we undertook 
that study, both analytically and numerically. Our point of departure was early theoretical 
work P-p!0|] on this problem, which we summarize for completeness in Sec. II. In particular, 
it had been predicted that a certain critical velocity exists at which the radiative power 



emitted by the kink would diverge. Below that critical velocity, radiation would be zero, 
and above it, it would decrease with increasing speed (see Ref. for a summary). Those 
results were obtained at a time where the main aim was to develop a perturbative approach 
to deal with soliton problems. That, and the fact that computers were not the easily 
available tool they are nowadays, meant that those results were never analyzed in depth 
or numerically checked. Therefore, as a first stage of our study, we devised a number of 
numerical experiments to check them, and we found no numerical evidence for the predicted 
divergence. In view of this result, we carried out an improved perturbative calculation, in 
the sense that it allowed us to interpret correctly the earlier results in Refs. [§-|TD| and to 
show that, although the earlier analyses were correct, the predicted divergence was actually 
unphysical. This theoretical analysis is reported in detail in Sec. Ill: A preliminary short 
report has been given elsewhere [0. The work done to that point suggested to us the idea 
that, opposite to what was beheved to date, sG kink dynamics on a periodic potential could 
be essentially that of a (pseudo-relativistic) particle. We thus applied a simple collective 
coordinate formalism to the problem, and it turned out to describe soliton behavior very 
accurately, even predicting unexpected new phenomena. The analytical approach and the 
numerical simulations are contained in Sec. IV. Finally, to complete our program studying 
length scale competition and its effects in nonlinear disordered systems, we performed further 
numerical experiments to clarify whether robust objects like kinks, which according to our 
collective coordinate theory behave mostly like particles, can still exhibit the destabilizing 
effect of length scale competition. We found that this was actually the case. Furthermore, 
the simplicity of kinks allowed us to carry out a (numerical) linear stability analysis which 
provided us with a clear explanation for the numerically observed features. We collect 
our results on this question in Sec. V. To conclude, we summarize the facts that we have 
learned, which considerably enhance the understanding of sine-Gordon kink propagation 
in disordered media and shed light in the so far unexplained phenomenon of length scale 
competition. Our results are also of relevance to many nonlinear systems of physical interest, 
mainly in three directions: First, all non- numerically validated or non-physically interpreted 
predictions obtained through perturbative calculations should be treated with a degree of 
caution. Second, the collective coordinate formalism yields a very simple and accurate way 
to deal with perturbed nonlinear problems, especially those in which the perturbation enters 
parametrically rather than additively. And third, length scale competition is an ubiquitous 
phenomenon that may be responsible for many instabilities arising in different nonlinear 
disordered systems. 

II. BRIEF SUMMARY OF PREVIOUS RESULTS. 

We start by describing the picture of sG kink propagation on parametric periodic po- 
tentials that has been accepted to date. The problem, that has been studied by Mkrtchyan 
and Shmidt p and Malomed and Tribelsky ||,|l3] is given by a perturbed sG equation of 
the form 

Utt — Uxx + [1 + e cos{kx)] sin u = 0, (1) 

(modeling for instance a long Josephson junction with modulated critical current, to mention 
just one application); the question posed was whether kinks can propagate freely in such a 



system, and if so, to describe this propagation. We will only record here a short summary 
of previous work. The reader is referred to the original papers ||^10| for details. 



Mkrtchyan and Shmidt ^ used a Green- function perturbation technique (GFPT). They 
derived a linearized equation for the first order correction to a kink moving with constant 
velocity, computed the Green function corresponding to its homogeneous version, and then 
used it to obtain the desired correction by integrating the source term with that Green 
function. They then noticed that radiation appeared only above a critical kink velocity 
^thr = (1 + k'^)~^^'^- At that particular value, the correction diverges, and the authors 
explain that their calculation become invalid in that region, as of course, it assumed the 
correction was small. On the other hand, the approach of Malomed and Tribelsky PJTO[] 
was quite different. Its basis was the Inverse Scattering Perturbation Theory (ISPT). A 
meaningful summary of this kind of calculation would be quite lengthy, and hence we will 
omit it here, referring the reader to Ref. p| which is mostly devoted to describing ISPT 
in detail. Let us just mention that the idea is that, if the amount of radiation emitted by 
the kink is small, as it should be if the perturbation is small, then the spectral density of 
the emitted energy can be computed following a Taylor expansion, and the total radiated 
energy is then derived by integration over all modes. Again, the result was that there was 
a critical velocity fthr = (1 + fc^)~^/^ such that kinks traveling with velocities v < Vthr 
did not emit any radiation at all, whereas in the opposite case the amount of emitted 
radiation decreased as f -^ 1, diverging when v = fthr- Always within the framework of 
ISPT, Malomed and Tribelsky |1^ were also able to show that dissipation could play a 



regularizing role, suppressing the divergence. As the results in Ref. agreed with those 
in Refs. [p|,p!0|, the existence of this threshold for radiation with its associated divergence 
was accepted, and the question of kink propagation on periodic potentials was regarded as 
basically solved. As mentioned above, it has to be borne in mind that the main issue of 
those early researches was to establish the proper foundations for a perturbative theory for 
solitons. Hence, the question of the physical meaning and origin of the divergence was not 
addressed. Another unexplained point arises already from ISPT, which allows computation 
of the radiation nature. When this is done in our case, the radiation wavenumbers turn 
out to be related to the perturbation one by a complicated equation (see, e.g., 0), which, 
in particular, implies that radiation is emitted with a non-intuitive wavenumber k~^ at the 
divergence. This prediction is difficult to understand physically. Let us recall at this point 
that a particlelike picture of kink propagation had been developed and had been largely 
successful so far [|I^ when compared to numerical experiments. If ISPT predictions for the 
radiation wavenumbers were true, the reason for them must come from the wave nature of 
kinks. Consequently, the particle picture should be regarded as a major simplification and 
valid only in limited cases. 

III. KINK PROPAGATION ON PERIODIC MEDIA. 

A. Numerical results. 

With the above scenario (and the question it poses) in mind, we carried out a num- 
ber of numerical simulations looking in the first place for the proposed threshold. All the 



simulations we will be reporting on throughout the paper have been carried out taking pe- 
riodic boundary conditions. The integration was performed with two different procedures, 
an adaptation of the energy conserving Strauss- Vazquez finite-difference scheme |T^ and 
a fifth order, adaptive stepsize, Runge-Kutta integration [Q of the discretization of the 



PDE. The results were independent of the procedure, which is a satisfactory checking. We 
performed a careful search, paying attention to the fact that the predicted value was a first 
order calculation, and that it may not be quantitatively accurate. On the other hand, the 
finite width of the simulated system may also be of relevance at this point, as its radiation 
spectrum structure is not identical to the continuum, infinite system (in particular, the low- 
est frequency in the model is restricted to be uj^^^ = 1 + {2tt/L)^, L being the length of 
the system). Hence, we monitored the amount of radiation emitted by the kink by making 
simulations with many different initial conditions, sweeping a range of initial velocities; if 
there was a threshold somewhere, there should be a change in the radiating power of the 
kink as it moved through it. The result was negative: No evidence for a threshold was found, 
even when the search was performed for a large range of initial velocities with a resolution of 
10^^ for some choices of k. Examples of the outcome of the simulations are shown in Fig. |1| 
for three values of the potential wavelength, of the order (a) and smaller (b, c) than the kink 
width (~ 6 in our dimensionless units) at f = fthr- It has to be stressed that the predicted 
divergence does not depend on the strength of the perturbing potential, e, but we also tried 
to make the effect more visible by increasing this parameter. Indeed, in Fig. |l], e = 0.4, a 
value that is not very small, and the kinks seem unaffected except for a small amount of 
radiation and an oscillatory motion superimposed on its trajectory, which is shown in Fig. 
0. It is interesting to note that the kink traveling on the short wavelength potential (c) ap- 
pears not to be affected at all. This will be understood by means of the collective coordinate 
approach in Sec. IV. On increasing e further, trapping behavior takes place, i.e., kinks are 
trapped by the potential and cannot propagate, but there is no strong emission of radiation 
(for an example, see Fig. ^(b), which will be discussed later). Actually, this trapping can 
be of two very different kinds, as we will discuss in sections IV and V. Another interesting 
remark is that we also observed that kinks always emit radiation, even when moving at a 
very low velocity, far below the predicted threshold. A similar result arises from the work of 



Peyrard and Kruskal on highly discrete sG systems ||15[, where kinks propagate on the peri- 
odic potential coming from the Peierls-Nabarro barrier, although this comparison should not 
be taken too literally as there are some differences between both problems, like the existence 
of a maximum allowed frequency in the discrete one, for instance. It thus becomes evident 
that the features of kink propagation on periodic potentials are qualitatively different from 
the above perturbative analytical results. Interestingly, numerical simulations on a similar 
perturbation of the 0'^ problem |jl6|] seem to confirm the absence of this divergence. We will 



elaborate more on this when presenting our conclusions in Sec. VI. 

B. Theory. 

In order to gain insight into the numerical observations, we developed a new perturbative 
approach for this problem, following a similar approach to that given by Fogel et al. |]T^ 



To this end, we perform a Lorentz transformation and rewrite (|l|) in the rest frame of the 
soliton (i.e., the reference frame moving with the speed of the unperturbed soliton, v) 



Utt — Uxx + {1 + ecos[A;7(x + vt)]} sin-u = 0, (2) 

with 7 = (1 — i;^)^^/^ the Lorentz factor. Here we consider e <C 1, so the perturbative term 
may be treated by assuming a solution of the form 

u{x,t) = Uv{x) + eu'^^\x,t)^ (3) 

where Uy{x) = 4tan~^(e^) is the unperturbed sG kink. For completeness, we now recall how 
the most appropriate basis in which to expand eu^^\x, t) is obtained. Introducing the Ansatz 
(§) in Eq. (^ without the perturbation term, linearizing in the small quantity u^^\x, t), and 
separating time and space by introducing u^^\x, t) = /(x)e~*'^*, we are left with the following 
eigenvalue problem for f{x): 



■-—- + fl - 2sech^x 



f{x)=cu'fix). (4) 



This is a well known eigenvalue problem |T^; there exists exactly one bound state, with 
Ub = 0, and a continuum of scattering states with u^ = I + k^; the corresponding normalized 
eigenfunctions are 

fb{x) = 2 sechx (5a) 

f{K,x) = ^e*"^ (/t + itanhx) . (5b) 



These eigenfunctions have a very clear physical meaning. The bound state fb{x) is associated 
to the Goldstone translation mode of the soliton, whereas the continuum eigenfunctions 
/(k, x) are the radiation modes (see [|I2| for a detailed discussion). Besides, these functions 



form an orthogonal basis, since the corresponding operator is self-adjoint. We will make use 
of this fact to deal with our problem. In terms of this basis, the first order correction can 
be split into two parts, namely, 

m(^) (x, t) = M^*''^'^^) (x, t) + u^'"^^ {x, t) , (6) 

where 

u(t--)(a;,t) = i0,(t)/,(x), (7a) 

o 

U^'^'^\x,t)= dK(P{K,t)f{K,x). (7b) 

J —oo 

To find the amplitudes 0fe(t) and 0(k, t), one again introduces the Ansatz (^ in Eq. (||), 
linearizes and Fourier transforms in time; subsequent projection yields 

(f)b{t) = 4 / dx cos[A;7(x + vt)] o— , (8a) 

j-oo cosh X 

, , 9x , / X f°° , r, / ,,e~*''^(K — itanhx) sinhx , ,, 

(f){K,t) + {l + K^)(f){K,t) = 2 dx cos[kj{x + vt)] ; 2—- (8b) 

J-oo J2tt(1 + k^) cosh X 



It now remains to solve Eqs. (P) and invert the various transforms needed to arrive at 
them. In the following, we discuss translation and radiative parts in (Bl) separately. Let us 



start with the simplest one, i.e., the translation mode contribution. Note that (|8a|) is, after 
performing the integration, nothing but the Newton's law for a time-dependent force. Its 
solution may be readily found, and finally one obtains 

«(t-s)(^^^) = ^ J sinik^vt) sechx. (9) 

Recalling that we are working in the unperturbed soliton reference frame, this is a localized 
oscillatory motion superimposed on its otherwise constant trajectory. Now, let us remark 
that the prefactor implies that short wavelength {k -^ oo) perturbations will have no effect on 
the motion of the center of the soliton, which is also in good agreement with our simulations 
in Fig. Ij. This behavior can be understood in terms of a "smoothing" of the potential: The 
kink, having a width much larger than the perturbation wavelength, experiences only an 
effective averaged force, whose amplitude vanishes exponentially for large k (see Sec. IV; see 
also related comments in 0,§]). 

Equation (|8E| ) for the K-mode radiative contribution can also be solved. After computing 
the integral in the right hand side of Eq. (^^, one is left with the Newton's law for a forced 
harmonic oscillator. This allows the determination of (f){K,,t), and substitution of it in Eq. 
( [TED to find the total radiative contribution: 



u^'-'\x,t)^- 



tanh X — -— 
ox 



(10) 



'■^ 1 + K^ - A;^7^ 



„ikjvt p—ikyvt 

+ 



cosh[7r(A;7 — k)/2] cosh[7r(fc7 + k)/2] 



^tnx_ 



It is possible to deal with the integral in (|T0]) in the complex plane: When x > 0, in the 
upper half plane, and when x < in the lower half plane. The pole structure of the integrand 
will completely determine the total radiative contribution. In particular, we will see that 
radiation only appears for some special values of the system parameters. 

We take x > in what follows (the opposite case is treated in the same way). Accordingly, 
the integral has to be analyzed in the upper half complex plane. The pole structure of the 
integrand is depicted in Fig. ^. All poles are simple, and their locations are zq = +i, 
Zi = ia = ^i\f\^^W^^y^ ^ and z^ = ±^7 + i{2n + 1), n being a non- negative integer. For 
the sake of clarity we treat each pole separately. 

i. The first pole, Zq = +i, is constant, and does not change when the system parameters 
change. Since this pole is purely imaginary, it is immediately seen that the contribution of 
the residue at Zq is exponentially localized around the kink center. This term does not give 
rise to any radiation, but rather to time-dependent corrections of the kink shape. 

ii. The family of poles z^ depends on the perturbation wavenumber k and on the kink 
velocity through the Lorentz factor 7. However, they always have a positive imaginary part, 
thus leading again to exponentially localized contributions. Therefore, the z^ poles also do 
not produce any true radiative correction. 

iii. The remaining pole is the key one. If a^ = 1 — fc^7^f^ > 0, the same reasoning 
applied to the other poles holds, and there is no radiation. It is worth mentioning that 
localized oscillations around the kink center, predicted from the contributions of zq, z\ [a 
real), and 2;^, were already evident in our numerical simulations, as shown in Fig. |l]. For 
fixed fc, as v increases, the pole moves down the imaginary axis, and at the critical value 



t>thr = (1 + fc^)~^/^ it lies at the origin of the complex plane. For kink velocities v > i^thr 
the pole is purely real, and then it does give rise to a radiative contribution, whose form is 
given by (with (3 = y/k'^'y'^v'^ — 1 a real number) 

(rad) TT 



U 



P 4^2^2 




i{k'yvt+l3x) —i{k^vt—(3x) 



cosh[7r(fc7 - /3)/2] cosh[7r(fc7 + /5)/2] 



:ii) 



This expression tells us that radiation occurs whenever /3 is real {v > fthr), and this radia- 
tion is the superposition of two linear waves of different amplitudes, travelling in opposite 
directions but with the same phase velocity. 

C. Discussion. 

To this point, it appears that our perturbative calculation leads exactly to the same 
prediction as those in |8|,p!0|, namely that there is a critical velocity Wthr = (l + A;^)"^''^ below 



which kinks do not radiate and above which they do. At that precise velocity, the amplitude 
of the emitted radiation diverges; notice that (3 vanishes as v approaches Wthr from above and 
consequently the prefactor in Eq. (|lT]) goes to infinity. However, this apparent equivalence 
is not so. The crucial difference arises when one looks more carefully at Eq. ([111) : As Wthr 
is approached, not only the amplitude of the emitted wave diverges, but also its wavelength 
2tt/j3. Then, we are faced with something similar to an "infrared" divergence, and usually 
those do not have a real physical meaning. We will show immediately that this is indeed 
the case here, but let us first comment on the reasons why our calculation provides us with 
this physically relevant result that was not transparent in the previous approaches. As to 
the GFPT computation [§, they compute the first order correction to the field much as we 
do here (actually the two approaches are basically the same in the beginning), but they do 
not use the natural translation mode-radiation basis, so they can not separate the different 
contributions and are therefore led to an expression they can not analyze in detail; as we 
already pointed out, they merely remark that their calculations are invalid in the vicinity 
of the divergence, as they assumed the correction should be small. On the other hand, 
ISPT ||^,|T^ yields a different result than ours in spite of using a suitable basis because the 
integration over k is made in an incoherent fashion, i.e., integrating over emitted energy 
instead of emitted amplitude (we notice in passing that many ISPT results are obtained by 
this same means). When the integration over radiation modes is made coherently as shown 
here, the result changes due to the superposition of different modes. These reasons lead us 
to believe that, although admittedly the early perturbative work was mathematically sound, 
the calculation we present here is the physically correct first order result. 

Now that we have a reliable perturbative calculation, we need to understand what is the 
nature of the divergence. To make progress, it is very important to turn to the form of our 
starting Eq. (|^) with dimensions, namely 

Utt — c^Uxx + i^o[l + ecos(A;a;)] sinu = 0, (12) 

where Cq and ujq are a velocity and a frequency characteristic of the particular physical 
context. Redoing the calculations with dimensions transforms the divergence condition 
^T^thr = 1 into fc7ofthr = ^0 [To = (l^^^/co)^"*^^^]- This immediately clarifies what happens: 

8 



The divergence occurs when the velocity of the kink is such that the time it takes to travel 
through a wavelength of the potential, Tq = A/(7fthr), A = 27r/fc, is exactly the period 
of the lowest frequency phonon, Tq = 2tt/ujq. If the velocity is lower than Vthr-, the kink 
will not be able to excite phonons, whereas when its velocity is higher it can and will 
subsequently radiate. At fthr, the excited radiation is that of the lowest phonon, and it has 
infinite wavelength and velocity, as predicted by our calculation. This natural picture of 
kinks exciting radiation according to the frequency of their propagation through a potential 
wavelength becomes therefore the likely candidate to explain the divergence. On the other 
hand, now it also becomes clear the divergence of the energy at fthr^ It diverges because 
of the infinite contribution arising from the infinite wavelength mode when integrated over 
the whole x axis. This agrees with GFPT and ISPT results whose only difficulty was not to 
specifically identify the mode responsible for the divergence. 

In spite of this clarification, the most significant question is not answered yet: Why 
numerical simulations do not agree with this calculation, which seems to allow for a simple 
and physically reasonable interpretation? By looking again at Figs. ^ and 0, it is easy to 
realize that the fiaw of the perturbative calculation is at its very root: We are computing 
first order corrections around a kink moving at a constant velocity v, and this condition 
never holds. Whatever the starting position of the kink is, it will behave like a particle in 
the sense that it will be accelerated or decelerated depending on whether it travels towards 
a minimum or a maximum of the potential. In fact, the translation mode correction itself is 
describing this: The kink velocity, in its reference frame, is not zero but rather it oscillates 
between positive and negative values. It is not a surprise, then, that first the resonance 
condition we have obtained is never matched, and second that the kink emits radiation at 
any velocity, because it is accelerating or decelerating. Of course, we should note that this is 
a perturbative calculation including only first order terms; the possibility still remains that 
the divergence is suppressed by higher order nonlinearities. 

IV. COLLECTIVE COORDINATE APPROACH. 

The above numerical results and the subsequent perturbative calculation strongly suggest 
that sG kinks behave as point-like particles in the presence of a periodic parametric potential 
like the one we deal with here. Therefore, it is natural to try to describe those results by 



means of the collective coordinate formalism. This approach was first proposed in |12], and 
it has been applied recently to sG breathers on periodic potentials [Q as well as to NLS 
equations with the same perturbation. In both cases the analytical predictions turned 
out to be in very good quantitative agreement with numerical simulations: For instance, 
in Ref. the threshold for breather breakup into a kink-antikink pair was predicted with 
an accuracy better than 0.1%. On the other hand, the calculation in Ref. predicted the 
appearance of the so called "soliton chaos," verified by simulations of the full PDE. In our 
present problem, the advantage we have is that, due to the simpler nature of the kink, we 
will be able to compute the effective potential not only for kinks at rest but also for moving 
kinks. 

The basic idea of the collective coordinate formalism is very simple: To reduce a compli- 
cated problem with an infinite number of degrees of freedom, posed in terms of a PDE, to 
a much less complex problem with a few degrees of freedom (and correspondingly described 



in terms of ordinary differential equations, ODE's). There are a number of ways to do this, 
and different quantities can be chosen to play the role of collective coordinates describing 
the motion of the nonlinear excitation as a whole. For our problem, it is enough to simply 
consider the kink center as our collective coordinate for the kink. Its motion will be then 
governed by an effective potential that can be computed by integrating the perturbative 
contribution to the hamiltonian over the kink profile, i.e., 

/oo 
dx[l — cosUv{x — Xo,t)]coskx, (13) 

-oo 

where u.u{x — xo,t) denotes now a kink moving with constant velocity v and centered at xq. 
This integral can be easily evaluated and yields 

Ves{xo, t) = 2e . — — - cos[fc(xo + vt)]. (14) 

7"^ smh(A;7r/27) 

From equation (0) we see that the potential experienced by the particle equivalent to the 
kink is basically the same perturbation potential that appears in the PDE ([l|), although the 
prefactor in front of it is quite complicated. The simplest dependence of this prefactor is 
on the wavenumber. It can be immediately seen that when k ^ (long wavelength limit) 
the effective potential prefactor reduces to 4/7 and subsequently Ves becomes closer to the 
perturbative one; in the opposite limit, k -^ 00, the sinh term makes the effective potential 
vanish exponentially. This is in agreement with what we have learned so far: Looking at Fig. 
|I], it can be seen that the short wavelength potential has no effect on the kink (c), whereas the 
motion on long wavelength perturbations resemble that of a particle on the bare potential. 
To phrase in the terminology introduced in Ref. [^], the behavior of the kink in these cases is 
that of a "bare" (long wavelength) or a "renormalized" (short wavelength) particle. It is also 
important to notice that this result agrees with the perturbative calculation we described 
in Sec. Illb [see Eq. (|)] as it was to be expected. There we showed that the correction 
to the center of mass motion was basically an oscillatory term, implying that the velocity 
of the center of mass oscillates around some mean value. This is precisely the same kind 
of trajectory followed by a point-like particle in the potential in Eq. ( |T4|) (at least if the 
velocity is not too close to 1). 

Nevertheless, it is worth pursuing this agreement a bit further, by studying the threshold 
for kinks to propagate in this kind of potential. The easiest way to compute the threshold 
is by equating the kinetic energy of the kink to the maximum of the effective potential, 
provided we restrict ourselves to the non-relativistic limit (f^ not too close to 1) to keep 
the kink mass constant. This will give us the maximum potential height over which a kink 
that starts from a point at which the perturbation is zero with a certain velocity is able to 
overcome the nearest top point. Using the fact that the mass of a not too fast kink is 8 in 
our units, we find that the threshold is given by 

ethr = —, smh --. (15) 

kiT 27 

In the same way, we could have computed the threshold velocity for a given strength of the 
potential, but we prefer to check our predictions this way because the presence of 7 makes the 
other possibility more complicated. We compared this prediction to numerical simulations. 
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We show an example of this comparison in Fig. H, where we study the propagation of 
a kink with initial velocity 0.2 (i.e., in a non-relativistic situation) on a potential of a 
moderately long wavelength. The predicted threshold for kink propagation as given by 
Eq. (|l^) is ethr = 0.0424; from Fig. H we see that the numerical result is bounded by 
0.043 < ethr < 0.0435, meaning that the error in our prediction is at most of 2% which is 
quite satisfactory. We have checked several other cases, which we summarize in Fig. |^; an 
excellent agreement is always found, even for very large values of the perturbation potential. 
To conclude this section on the collective coordinate treatment of the problem, we discuss 
another prediction of it that it is numerically verified, relativistic effects playing now the 
relevant part. By looking again at Eq. (|I4|), it can be realized that the presence of 7 in 
the potential may give rise to singularities in the neighborhood of the maximum velocity, 
V = 1. To check whether this is so, we integrated numerically the ODE obtained from 
applying second Newton's law to Ves, and found that if the initial conditions were those 
of a particle starting at the top of a large potential, the velocity of that particle would 
grow as it slides down the potential; of course, if the potential is large enough, the velocity 
can reach v = 1: In those cases the numerical integration of the ODE broke down. The 
question then arises whether this is an artifact of our collective coordinate approach or there 
is a related phenomenon in the full PDE. We actually found that this ODE prediction is 
verified, as we shown in Fig. ^ In this simulation, the initial condition was a kink at a 
maximum of the potential with initial velocity 0.1 so as to start the motion from the stable 
point. As it moves to the nearest well, it accelerates and, eventually, its velocity becomes 
very close to 1, implying that the kink cannot accelerate further. Then it is trapped at an 
intermediate point of the potential, instead of continuing its motion to the bottom. The 
existence of this counterintuitive phenomenon shows in a very dramatic way the value of 
a simple approach like the collective coordinate formalism to help understand complicated 
nonlinear phenomena. 

V. LENGTH SCALE COMPETITION. 

A. Numerical experiments. 

The numerical and theoretical analysis described so far allowed us to achieve a quite good 
understanding of the periodically perturbed sG problem, at least of the basic phenomenology. 
With that background in mind, we then turned to the initial motivation for this work, 
namely to study kink propagation in periodic potential as a step towards the much more 
complicated problem of breather propagation on periodic media ||^. In principle, we did not 
expect length scale competition to arise in this problem, as kinks are very robust objects 
(which is further confirmed by our above results, in particular by the success of the collective 
coordinate approach) . On the other hand, opposite to the case of the breather, kink widths 
do not vary much when changing the only parameter governing it, the kink velocity. Of 
course, when approaching the maximum velocity, Lorentz contraction of the kinks will make 
them vanishingly small, but that is a regime in which is very difficult and time consuming 
to carry out good numerical simulations, so we did not address the problem in that limit. 
Therefore, the kinks we are usually dealing with have a width of about 6 in our dimensionless 
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units. Our purpose was to perform some simulations in the high perturbation regime to see 
whether any hght could be shed on the related breather problem. 

The numerical experiments we made were as follows. We studied several potentials of 
different wavelengths, ranging from 0.5 to 20, i.e., from much smaller than the kink width 
to roughly three times its width. The initial condition was always a kink at one top of the 
potential; different velocities were considered. A summary of these experiments in the most 
interesting range of potential wavelengths is shown in Figs. |^ and ^. Figure ^ shows the 
motion of the center of mass of the kink on different potentials. The motion in the small 
or large wavelength limits has been already discussed and it is again seen here. However, 
a more interesting phenomenon is also evident, namely kink trapping [or even reflection in 
the case of wavelength 3, see Fig. H(b)]; this trapping was not to be expected because kinks 
start from a maximum of the potential and with a large initial velocity (in Fig. |^, it was 0.5). 
Regarding this point, we have to stress here that this trapping is of a different kind that the 
one discussed in the preceding section, which was clearly a non-radiative process. Besides, 
the trapping depends crucially on the wavelength of the potential. Thus, for instance, in the 
case of wavelength 2 [Fig. §|(a)], the kink is able to propagate over six wells, whereas when 
the wavelength was 4 or 5 [Figs. |Kc), (d)], it was trapped already on the second well. This 
dependence clearly indicates that length scale competition is to same extent present also 
in the kink problem. This hint is further supported by our previous work for the breather 
case [0, which showed that this competition was most effective when the wavelength of the 
perturbation was around half the breather width or slightly larger. This is also the case in 
these simulations. Another common feature between both problems is that the outcome of 
a simulation depends very sensitively on the initial condition. This can be understood from 
the reflection case in Fig. ^(b): For the kink to jump back over one potential wavelength 
it is necessary that it meets the radiation it left behind in the appropriate phase to gain 
energy from it, and this evidently depends crucially on the initial velocity, as we checked 
in our simulations. Hence, we conclude that these phenomena are a novel manifestation of 
length scale competition. 



B. Linear stability analysis and discussion. 

The numerical findings we have described in the preceding subsection are of great impor- 
tance: The existence of processes governed by length scale competition in the sG kink case 
opens the possibility of a deeper study of the mechanisms through which this competition 
affects the kink evolution. The relevance of this comes from the fact that, when studying 
the sG breather problem [^], we were not able to carry out this deeper analysis due to the 
more complicated nature of the breathers, namely their intrinsic internal dynamics which 
severely complicate linear stability analysis. But after showing that this competition also 
affects kinks, we can certainly study their linear stability analysis, and, consequently, obtain 
for the first time an understanding of the mechanisms underlying length scale competition. 

We tested the stability of the analytical continuum solutions as well as of the numerical 
solutions in the following way. Let the solution to discretization of the sG equation (|I]) be 
Ui = Ui + Vi, where u\ is either the discretized version of the continuum kink or the true 
minimum energy static solution of the perturbed sG equation (^ obtained numerically, and 
where Vi is a small discrete- valued function whose time dependence is given by sin(u;t); the 
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index i runs over the A^ points of the discrete lattice. The discrete version of the perturbed 
problem (|l|) is given by 

Ui — a^^('Uj+i — 2ui + Ui-i) + [1 + ecos(A;ai)] Ui = 0, (16) 

where a is the lattice spacing. Substituting the proposed form for Ui in Eq. (|16D and lin- 
earizing we get 

nv = to'^v, (17) 

where v is the vector containing the Vi, i = 1, . . . ,N and f2 is an A^ x A^ matrix given by 

{2 + [1 + e cos(A;ai)] cos (u\ j , if i = j, 
-1, if^=J±l, (18) 

0, otherwise. 

In this formulation, the modes for the linear excitations around the shape u^^^ are obtained 
simply by solving for the eigenvalues a;^ of the matrix f2. We did this for all the wavelengths 
we were studying, taking for u^^^ the exact continuum sG kink at the top or at the bottom 
of the potential, as well as the numerically obtained solutions at similar positions. 

The results for our numerical linear stability analysis are shown in Fig. |^ for some 
of the relevant wavelengths |l^ . There are a number of interesting features which deserve 
comment. First, let us consider the spectra for the exact continuum shapes. When placed at 
the top of the potential, this gives rise to a negative lowest eigenvalue u"^, indicating that this 
continuum kink is not an exact solution of the discrete problem and has a tendency to relax to 
the correct one by emitting a burst of radiation. The shape at the bottom of the potential 
does not show this negative eigenvalue but, instead, a single bound state with frequency 
uj"^ ~ 0.3. This corresponds to a shape mode, similar to that present in unperturbed (j)'^ 
kinks, and it actually shows up in simulations: Initial data placed at a potential well (an 
exact continuum kink) exhibit a static center of mass but a general oscillating shape. This 
is easily understood if one realizes that in this range of potential wavelengths, different 
parts of the kink undergo the action of very different perturbation values (that can be even 
positive or negative contributions). In response to these gradients, the continuum kink 
oscillates. These isolated states characteristic of the continuum kinks disappear when we 
analyze static numerical solutions. It is seen from Fig. |^ that in that situation spectra are 
composed of bands. Actually, this is a common feature to all analyzed shapes, including the 
continuum ones, and it could be expected in view of the following argument: Far from the 
kink center, which only spans a small fraction of the lattice sites, the nonlinear contribution 
to the linearized discrete problem for Vi vanishes, and one is left with what is a standard 
Floquet (Bloch) problem. The corresponding structure is very well known, and in fact it is 
very much like the ones we show here, with gaps at positions that depend on the potential 
wavenumber and gap amplitudes that depend on the potential strength (e). 

By comparing the spectral structure we have obtained to the outcome of the numerical 
simulations (Fig. |^) and especially to the radiation emission (Fig. |^), the mechanism for kink 
destabilization can be inferred. In the case of small wavelength potentials, the spectrum is 
very similar to the unperturbed sG one (see Fig. ^). Therefore, the behavior of the kink is 
very similar to the continuum one moving in a discrete lattice, the periodic potential then 
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being nothing but the Peierls-Nabarro barrier, as we aheady mentioned. In that case the 
kink is known to radiate [^ and correspondingly decelerate until it is eventually trapped 
in a potential well (Fig. ^d). When the wavelength is smaller, the height of the effective 
potential seen by the kink is so small [due to the sinh term in Eq. (|14D] that this effect is 
hardly noticeable (hence the perfect constant motion of the kink on the smallest wavelength 
potential in Fig. |^). As the wavelength increases, more and more modes move below the 
phonon band, inducing shape oscillations of the kink, and as a consequence of this motion, 
long wavelength radiation is emitted [clearly seen in Figs, ^(a) and (b)]. This is possible 
because in those cases there is still a large number of available modes just above the phonon 
band. Note that the lower limit of the phonon band is given by cu = 1; lower frequencies 
are localized, because they cannot propagate in the system far from the kink, where the 
spectrum structure is essentially the unperturbed one. This combined effect induces a rapid 
destabilization of the kink and its trapping. Finally, if the wavelength is further increased, all 
the first band eventually moves below the phonon band (Fig. ^d) , and hence long- wavelength 
emission is strongly suppressed (Fig. ||d), which stabilizes the kink, making possible its 
propagation. The effect of the shape modes coming from the first band is still revealed by 
the kink shape oscillations [see Fig. P(d)]. We believe that this interpreation clearly explains 
the mechanisms leading to the appearance of length scale competition in sG kinks. To seek 
further evidence, we looked for the approximate value of the potential wavelength at which 



the last mode in the first band crosses the phonon band. It is shown in Fig. |10| that this 
happens for a potential wavelength between 5.2 and 5.3. Were our hypothesis true, kinks 
would not be able to propagate on the former potential and they would be able to do it 
in the latter one. The numerical simulations shown in Fig. |Tl] confirm that this is indeed 
the case. Interestingly, the radiation is quite different in both cases and, furthermore, for 
the trapped kink trapping occurs at the second potential well, indicating that the small 
number of modes available to radiate prevents a very rapid decay of the kink. We thus 
conclude that our interpretation is indeed correct. As a matter of fact, as this feature of 
the spectrum will also be present when considering breather propagation (recall that the 
reason for the appearance of the gaps is the perturbating potential acting on the wings 
of the excitation, and the nature of the center becomes less relevant), our explanation of 
length scale competition should also apply to breathers. The results in Ref. |Q are in perfect 
agreement with what we have described in this section. 

VI. CONCLUSIONS. 

In this paper, we have studied kink propagation in ID sG systems with a spatially peri- 
odic modulation of its characteristic frequency. We have shown numerically that kinks can 
propagate steadily and mostly undisturbed, even for large amplitudes of the perturbation. 
Disagreement with analytical predictions previous to this work is resolved through a new 
perturbative calculation. By this means, the radiative divergence is shown to be similar to 
an "infrared" divergence. Our calculation provides us with a good physical understanding 
of the problem of free kink propagation in the periodically modulated sG system, which had 
not been obtained from the previous ones, of a more formal character. A comment is in order 
here, regarding the fact that now the mathematical foundations of perturbation theory for 
solitons are already established (mostly by pioneering works like Refs. [§,0), the emphasis of 
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that perturbation theory should be focused on their physical implications. Therefore, per- 
turbative calculations in nonlinear equations should be regarded as speculative if they are 
not verified through numerical simulations, and, most importantly, if their physical mean- 
ing is not fully established. That is one of the most important points of this work: After 
an appropriate perturbative calculation, and by carefully considering the dimensions of the 
problem, we have been able to identify the underlying physical reason for that divergence 
as a resonance with the lowest, infinite wavelength phonon mode. Moreover, motivated by 
our perturbative results, we have developed a collective coordinate approach to this problem 
that describes in a quantitatively correct way the main features of kink propagation, show- 
ing that the already known ^ "bare" and "renormalized" particle limits apply also in this 
case. The collective coordinate equations turn out to predict counterintuitive phenomena 
whose existence is confirmed by numerical simulations, namely kink trapping at intermedi- 
ate points in the potential. This is an important success of the technique. Finally, we have 
shown that length scale competition arises unexpectedly in this problem, which afforded us 
the opportunity to increase our understanding of this ubiquitous phenomenon [^|[ . By a 
detailed linear stability analysis, we have identified the mechanism by which length scale 
competition arises as coming from the band structure induced by the perturbation poten- 
tial. Again, the predictions of our theory have been fully confirmed (quantitatively) by the 
corresponding numerical simulations. 

The global picture that emerges from this work is that, once again, sG kinks behave 
basically like particles and a collective coordinate approach can be more faithful than com- 
plicated perturbative results. Length scale competition phenomena are only relevant in a 
large amplitude perturbation regime and therefore do not interfere with this picture in most 
situations. In this respect, it has to be noted that a perturbative calculation describes ev- 
erything beyond the center-of-mass dynamics: Extended (background) contributions (see, 
e.g., the third reference in |jl2[); shape changes localized around the (moving) kink; and the 
radiation, i.e., emission from the kink. It is crucial to separate and identify these physically 
different effects if one is to properly understand the dynamics of the considered system. Our 



results are likely to be general for kink-bearing models in view of the related results of |T^ 
on the spatially periodically perturbed 0^ model. Some remarks are in order regarding this 
related problem. The same kind of divergence is predicted by a pertubation theory similar 



to the one used here (see |T^ for details on this approach), and again, numerical simulations 



show evidence of the unphysical character of the divergence: It can be seen from Fig. 4 of |T^ 
that as the velocity of a decelerating kink goes through the threshold nothing special occurs. 
As a matter of fact, most of the reasoning we have used in the present study applies also 
to that work (with an additional feature coming from the shape mode of the unperturbed 
0^ kink), thus reinforcing the generality of our results. On the other hand, the discovery of 
length scale competition for sG kinks has allowed us to understand the underlying physical 
reasons, and helped us to gain insight on related results for the sG breather [Q . We believe 
that the mechanism we have identified in this work is of a very general character, and it 
will be important to have more work on related systems to check our predictions. As a final 
remark, we want to stress that we have provided a quite thorough description of the fea- 
tures of free kink propagation in a periodically modulated sine-Gordon system, and that we 
have been able to provide a consistent, physical framework to understand this "canonical" 
problem, that will be of help in dealing with more complex situations. 
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FIGURES 

FIG. 1. Absence of radiative divergence for kinks propagating in the spatially periodic sG 
model. Parameters are: (a) k = 27r/5, initial velocity v = Wthr = 0.387726. . .; (b), k = tt, initial 
velocity v = vthr = 0.091999 . . .; (c), k = 27r, initial velocity v = t;thr = 0.024704 . . . [corresponding 
wavelengths are (a) 1, (b) 2, (c) 5]. In all three cases, e = 0.4. The amplitude of the emitted 
radiation is very small; due to the periodic boundary conditions, it can be seen reentering the 
simulation interval without any appreciable interaction with the kink. Only half of the simulation 
interval is shown in plots (b) and (c) to enlarge details. Time increases upwards with final time 
t = 100. The potential is indicated by the dashed line (amplitude not to scale). 

FIG. 2. Instantaneous center of mass positions as obtained from the simulations in Fig. |^. 
Dot-dashed line: k = 2it/5. Dashed line: k = ir. Solid line: k = 2tt. 

FIG. 3. The pole structure of the radiation contribution. Filled circles mark the location of 
the poles which give rise to corrections localized around the soliton. Empty circles denote the 
locations of pole zi as the velocity changes. Only when zi becomes real {v > fthr) does it originate 
propagating wavelike corrections. See text for further explanation. 

FIG. 4. An example of the verification of the collective coordinate approach predictions. For 
a kink starting from a mid-point of a potential of wavelength 20, with velocity 0.2, the predicted 
threshold for propagation is ethr = 0.0424. (a) e = 0.43 and the kink propagates; (b) e = 0.435 and 
the kink is refiected by the potential maximum; (c) center of mass motion for better comparison 
of both cases; solid line corresponds to the simulation in (a) and dashed line to that in (b). Final 
time is t = 200. Notice the absolute absence of radiation in this phenomenon. 

FIG. 5. Comparison of the effective potential prediction with numerical simulations for a num- 
ber of cases. Solid line: Threshold for initial velocity 0.2 as a function of the potential wavenumber. 
Dashed line: The same for initial velocity 0.5. Diamonds and crosses are points obtained from the 
numerical study. Error bars in the numerical threshold are of the order of the size of the symbols. 

FIG. 6. Trapping of a kink at an intermediate position in the potential instead of at the bottom. 
e = 2 and the kink starts with velocity 0.1. The amount of emitted radiation is large due to the 
fast acceleration of the kink in this large perturbation, (a) time evolution of the kink; (b) time 
evolution of the center of mass. The neighboring minimum is at x = —2.0 (indicated by the dashed 
line), and the kink oscillates around approximately x = —2.5. Final time is t = 50. 

FIG. 7. Length scale competition. Center of mass evolution of a kink on potentials of different 
wavelengths, always with e = 0.7 and initial velocity 0.5. Wavelengths are: O, 1; -|-, 2; □, 3; x, 4; 
A, 5; and *, 6. See text for explanation. 



FIG. 8. Time evolution for some of the kinks in Fig. |^ Wavelengths are (a) 2; (b) 3; (c) 4; and 
(d) 6. Notice the different kind of radiation emitted in each case. Notice also that the kink in (b) 
goes back over one maximum of the potential probably due to resonant interaction with radiation 
left behind. Final time is i = 50. 

FIG. 9. Frequency spectra for kinks on periodic potentials. See text for details on how is it 
computed. Spectra for numerical shapes are plotted with O (potential maximum), and □ (potential 
minimum). Spectra for analytical shapes are plotted with + (potential maximum) and x (potential 
minimum). Wavelengths are (a) 2; (b) 3; (c) 4; and (d) 6. Notice in this last case that all the first 
branch of spectrum is below the phonon band. 

FIG. 10. Spectra for kink on periodic potentials of wavelength 5.2 (O) and 5.3 (x). Between 
this two values all the first band of spectrum occurs below the phonon band. 

FIG. 11. Numerical verification of the threshold for length scale competition. Evolution of a 
kink on a potential of wavelength (a) 5.2, or (b) 5.3. (c) shows the motion of the center of mass in 
both cases for better comparison; solid line corresponds to the simulation in (a) and dashed line 
to that in (b) 
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